Field-Induced Transition on Triangular Plane in the Spin Ice Compound Dy 2 Ti 2 07 
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The origin of the lowest-temperature anomaly reported several years ago using a polycrystalline 
sample of the spin ice compound Dy2Ti2C>7 had remained unresolved. Here we finally clarify its 
origin by susceptibility measurements down to 65 mK using single crystals under accurate control 
of the magnetic fields in two independent directions. We demonstrate that the transition is induced 
under a subtle field combination that precisely cancel the nearest-neighbor spin interactions acting 
on the spins on the triangular lattice within the pyrochlore structure. Contrary to the other two field- 
induced transitions, this transition is driven only by the interactions beyond the nearest neighbors. 
Our observation thus provides the first qualitative evidence for the essential importance of the 
dipolar interaction beyond the nearest neighbors in the spin ice. 

PACS numbers: 75.30.Kz, 75.40.Gb, 75.50.Lk 



Recently, the physics of geometrically frustrated 
systems has attracted a great deal of attention, because 
a variety of unconventional spin behavior emerges as a 
result of the suppression of ordinary long-range ordered 
states. Among them, the spin ice behavior is peculiar 
since the spin frustration originates from ferromagnetic 
(FM) interactions, in contrast with frustration due to an- 
tifcrromagnctic (AFM) interactions commonly found in 
systems based on triangular spin configuration. In spin 
ice compounds such as Ho2Ti 2 07 and Dy 2 Ti 2 07 0, 
crystal field effects endow the ground state of Ho 3+ and 
Dy 3+ ions with a doublet with strong Ising anisotropy 
along the local (111) direction Q. Since these ions have 
large magnetic moments, the dominant interaction is FM 
dipolar interaction D n n, not AFM exchange interaction 
J nn , and the effective nearest neighbor (nn) interaction 
Jeff (= flnn + Jnn) is FM. Throughout this paper, we 
use the term "spin" to designate the magnetic moment 
of a rare earth ion. The FM-nn interaction leads to 
the "2-in 2-out" ground state in which two of the four 
spins on each tetrahedron point inward and the other 
two outward. Such a "2-in 2-out" constraint, called 
the ice rule from the analogy to the proton ordering in 
water ice, and the resulting spin configurations induce a 
macroscopically degenerate ground state with a residual 
entropy. 

At first, the spin ice behavior was expressed by the 
near neighbor spin-ice model considering only J g 
However, in order to explain some key experimental 
features quantitatively the dipolar spin-ice model was in- 
troduced, in which the importance of long-range dipolar 
interaction is emphasized Q. By Monte Carlo simula- 
tions based on this model, it was possible to reproduce 
the results of powder neutron diffraction patterns 0. 
The model further predicts a phase transition to the 



true ground state using a spin loop move algorithm |6j, 
although the predicted transition has not been observed 
experimentally down to 65 mK Q. In addition, it has 
recently been shown that long-range dipolar interactions 
are largely screened but the remaining medium-range 
interaction within several neighbors plays an important 
role in the (/-dependent susceptibility and is the cause 
of the ordering ||. For definitive evidence of the 
importance of the dipolar interaction apart from the nn 
interaction, it is crucial to demonstrate a phenomenon 
induced directly by the further neighbor interactions. 
An outstanding problem related to the above discussion 
concerns the origin of the peak in the specific heat 
divided by temperature C/T around 0.35 K observed 
in polycrystalline Dy 2 Ti 2 07 by Ramirez et al. 0]. 
They reported three peaks at 1.2, 0.45 and 0.35 K in 
C /T. Curiously, the peak temperatures are almost 
independent of the applied magnetic field strength. 
From the specific heat for the [110] field direction |j, 
the origin of the peak at 1.2 K has been ascribed to the 
ordering of spin chains with the Ising axes perpendicular 
to the field direction [ToL HlT |. In contrast, the origin 
of the 0.45 K peak is most probably related to the 
transition from "Kagome-ice state" to 1-in 3-out ordered 
state induced when magnetic field is parallel to the [111] 
direction 0, 0, Q, 0|. However, to the best of our 
knowledge there have been no experimental reports 
clarifying the origin of the lowest temperature peak at 
0.35 K. 

Very recently, J. Ruff et al. theoretically proposed 
an origin for the 0.35 K peak 01- The pyrochlore 
lattice can be viewed as an alternating stacking of 
the triangular and Kagome planes normal to the [111] 
direction (Fig. 1(a)). In a range of field direction within 
the Kagome plane centered at the [112] direction the 
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FIG. 1: (color online) (a) The relation between the pyrochlore 
lattice and two field directions. A pyrochlore lattice consists 
of an alternating stacking of the triangular (red) and Kagome 
planes (blue) along the [111] direction. The [112] direction is 
in the Kagome plane, (b) One of the possible spin configura- 
tions of a pair of tetrahedra. The arrows at the vertices of a 
tetrahedron represent the magnetic moment of Dy 3+ ions. 



directions of all the spins on the Kagome planes are 
determined by this field. However, the spins on the 
triangular planes are perpendicular to this field and 
are experiencing no Zeeman interaction. Thus, the 
direction of these spins is decided only by the effective 
nn interaction, J G ff. As shown in Fig. 1 (b), three 
spins (1, 2 and 3) on the Kagome plane are fixed in 
a 1-in 2-out configuration for the up-tctrahcdron and 
2-in 1-out for the down-tetrahedron. Then the ice-rule 
constraint makes all the spins on the triangular planes 
(4) point along the [111] direction. In this situation, 
an additional [111] magnetic field may cancel out J c ^ 
on spins labeled (4). These spins (4) are among them- 
selves third nearest neighbors on the pyrochlore lattice. 
Combined spins on different triangular layers form a 
face centered cubic (FCC) lattice when effects from the 
spins on the Kagome planes are neglected. Under such 
a canceling field, the residual spin-spin interactions, 
which mainly consists of long-range dipolar interaction 
and the third nearest-neighbor exchange interaction J3, 
become important. Such interactions are expected to 
induce a FM ordering of these spins if J3 is not too 
strongly AFM [ll|. Here we experimentally identify the 
predicted field-induced transition as the origin of the 
anomaly at low temperature. 

We used single crystals of Dy 2 Ti 2 07 grown by a floating 
zone method Q, and subsequently annealed at 700 °C 
in O2 for 24 hours. Two samples were cut from single 
crystalline rods into a thin plate with the plane normal 
to the [111] direction and a pencil-shape rod parallel to 
the [111] direction. The size and mass of the two samples 



were approximately 1.7 x 1.7 x 0.7 mm 3 and 14 mg, and 
0.8 x 1.0 x 3.2 mm 3 and 19 mg, respectively The [112] 
direction lies on the (111) plane. The real component 
X' and the imaginary component x" °f ac susceptibility 
were measured by a mutual inductance method down 
to 65 mK using a commercial dilution refrigerator. The 
strength of the ac field was 0.05 Oe-rms along the [111] 
direction. In order to control the magnetic field along 
the [112] and [111] directions independently, we used a 
"vector magnet" system consisting of two superconduct- 
ing magnets 0. One of the magnets points along the 
horizontal axis H x and the other one along the vertical 
axis H z (aligned to the [111] direction). We can apply 
fields up to 5 T for H x and 3 T for H z . These two mag- 
nets are rotated around the cryostat using a pulse motor 
with an angular precision of 0.01 degree. To ensure good 
thermal contact, we evaporated gold films on the sample 
surfaces and linked the sample to the heat bath using 
a 100-/im gold wire. We used a pair of Hall sensors for 
obtaining an accurate field alignment along the [112] 
and [111] directions. We determine the azimuthal angle 
0k in the Kagome plane between the directions of H x 
and [112] by the signal of Hall sensors; at </>k = 0, H x 
is parallel to the [112] direction. The overall accuracy 
of the field alignment including the resolution of the 
Hall sensors and misalignment between the Hall sensors 
and the crystal axes were within 2 degrees for the [112] 
direction and 1 degree for the [111] direction. Because 
the field direction changes continuously, it is difficult 
to estimate an accurate demagnetization factor, and we 
present data here with no demagnetization correction. 

In the measurements, we first applied fioH x = 
3 T along the [112] direction in order to hold the 
spins on the Kagome planes and then varied an ad- 
ditional H z along the [111] direction. At hqH x = 
3 T, the Zeeman energy for the spins 1 (Fig. 1 (b)) is 
£?2 = gjJ^{^H x )cos6i = 19.1 K and for the spins 2 

and 3, = Ef> = gjJ/j, B (Mo#*)cos0 2 = 9.6 K. Here 
gj is Lande's g factor, J = 15/2 is the total angular 
momentum, fi B is the Bohr magneton, fi is permittivity 
and 9i is the angle between the field direction and 
Ising axis at site i. These values are much greater 
than both the effective nn interaction, J e g = 1.11 K £| 
and the measurement temperatures in this experiment. 
Therefore, at relatively small H z /H x , all the spins on 
the Kagome planes arc fixed by H x . 
In Fig. 2, we show the temperature dependence of 
x' and x" & t l^oH x = 3 T, at various /-IqH z and at 
497 Hz for <f>K = 0° for a pencil-shape sample. Both 
components are highly asymmetric with respect to H z . 
At n$Hz = 0.20 T, the changes are most pronounced. In 
contrast, at fioH z = - 0.20 T and at other values of H z 
both components change little with temperature. Such 
enhancements around 0.2 T are also observed clearly 
in H z dependence discussed later (Fig. 3). In addition, 
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FIG. 2: (color online) Temperature dependence of real and 
imaginary parts of the ac susceptibility (a) x' an< i (b) x" 
of Dy2Ti207. Note that both components are highly asym- 
metric in H z - Features below 0.3 K for fioH z = 0.20 are 
attributable to a magnetic transition. 



even at higher temperature (~ 0.8 K), both components 
are enhanced at noH z = 0.2 T. This large response to 
the small ac field strongly suggests the spins on the 
triangular planes are freed from the spin-ice constraint 
and the external fields, and become paramagnetic. This 
is exactly expected by the theory of Ruff et al. 01 that 
the small [111] field (~ 0.24 T) cancel out J e g for spins 
on the triangular planes. 

At [IqH z = 0.20 T, there is a broad peak centered at 
0.29 K in \' and a peak at 0.20 K in \" ■ We also 
observed these peaks at almost the same temperature 
T = 0.28 ± 0.01 K for different frequencies (50, 1.75 k 
and 3 kHz; data not shown). This invariance of the peak 
temperature indicates that the peaks are associated with 
a transition to long-range magnetic order rather than 
with spin freezing |l7| ■ The results of thin plate sample 
are qualitatively the same as those of pencil-shape 
sample (data not shown). 

In Fig. 3 we show the H z dependence of \' at noH x = 4 
T and at 0.30 K for K = 0° and 180°. Most importantly, 
we find that at 0k = 60 x2n deg. (n = integer), \' 
exhibits a peak at a certain positive H z (fioH z ~ 0.2 T, 
0k = 0°). In contrast, at 0k = 60(2n+ 1) deg. the peak 
in x' appears at the corresponding negative H z (pqH z 
~ - 0.2 T, K = 180°). This field reversal is indeed ex- 
pected by the difference of the stable spin configuration 
under different H x direction. When we rotate H x by 60 
degrees, the stable spin configuration on the Kagomc 
lattice changes from 2-in 1-out to 1-in 2-out for the 
up-tetrahcdron and vice versa for the down-tetrahedron 
(Fig. 1(b)). At the same time, the stable spin direction 
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FIG. 3: (color online) Dependence of the real part of the ac 
susceptibility \' on the field strength along the [111] direction 
for (f> K = 0° and 180°. When H x is rotated by 180 (60 x 3) 
degrees, the sign of the peak field changes. 



on the triangular planes changes from the [111] to [111] 
direction and the sign of the required canceling field H z 
changes. The observed asymmetry in H z and the sign 
reversal of the characteristic H z by the rotation of H x 
strongly suggest that this anomaly is indeed due to the 
predicted spin ordering on the triangular planes 
If the field alignment is accurate, the absolute value of 
the peak field H z should be the same. In reality, there 
is a little misalignment of both fields in the experiment. 
A small difference among the peak fields at different 
0k is attributable to a small misalignment of H x from 
the (111) plane. From the H x dependence of the peak 
field (data not shown), we can evaluate the critical field 
by extrapolating the peak field in H z for various H x to 
fJ'oHx = T and the misalignment angle by considering 
the difference between the critical field and actual peak 
value. From this procedure using the data at n = to 
5 for 0k, we evaluated the critical field fioH c as 0.25 ± 
0.01 T and the misalignment as about 2 degrees. This 
H c is consistent with the theoretical expectation. The 
expected Zeeman energy of a spin on the triangular 
planes by the combined nn exchange and dipolar energies 
is 1.71 K considcringlong-rangc dipolar interaction by 
the Ewald method [llj. Thus, the scale of the critical 
field H c is estimated as 1.71 K/{gj{J z )ii B /k B ) = 0.25 T. 

In Fig. 4, we show the H z dependence of x' and x" at 
0.3 K and at various H x for 0k = 0°. With decreasing 
horizontal field H x , the peak around /ioH z = 0.22 T 
disappears below [ioH x = 1.2 T. This is compatible with 
the behavior of C/T of polycrystals reported by Ramirez 
et al. 0] for which the peak appears only for fields above 
1 T. This provides additional evidence that the origin of 
the low temperature peak in C/T is the ordering of the 
spins on the triangular planes. 
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FIG. 4: (color online) Variation of the ac susceptibility peak 
with the field H z along the [111] direction and H x along the 
[112]. The peak around noH z — 0.22 T disappears below 
HoH x = 1.2 T. 



These anomalies in x' and x" emerge under the condition 
that the effective nn interaction for the triangular spins 
are canceled out and the long-range dipolar and further 
neighbor exchange interaction between these spins are 
important. Until now, there have been some discussions 
about the importance of the long-range interactions for 
explaining the physics of spin ice, for example, fitting 
the specific heat C v 0, and the neutron structure 
factor S(q) 0. In these discussions while of quantitative 
relevance, the importance of the long-range interaction 
was not made so qualitatively clear and important. Our 
report here marks the first experimental demonstration 
of a phenomenon in the spin ice originating from the 
long-range interaction. 

The observed critical field is in good agreement with 
the theoretical expectation. The predicted process also 
explains the behavior in polycrystals: grains satisfying 
the critical condition /Zo-ff[ui] = Mo-^c = 0.25 T with 
/io-ff[n2] > 1-2 T exhibit magnetic transitions one after 
another with increasing magnetic field. In the Monte 
Carlo simulation using a simple dipolar spin ice 

model with long-range dipolar interaction and only the 
nn exchange interaction J nn , FM ordering of the spins 
on the triangular planes occurs at 0.65 K. When the 
AFM third nearest-neighbor exchange interaction J3 
= - 0.027 K is included (J3 = in the simple dipolar 
spin ice model), the transition shifts to 0.28 K with 
the ordering of the spins on the triangular planes still 
FM. Although it is difficult to conclude solely from the 
present ac susceptibility measurement that this ordering 
is FM ordering as theoretically predicted, the observed 



large peaks in the dissipation x" are consistent with the 
motions of FM domains. 



The estimated value of J 3 from this procedure, how- 
ever, is not consistent with that to reproduce the transi- 
tion under the [110] field, which is well fitted with J 3 = 0. 
To resolve this discrepancy, one may have to reconsider 
the exchange parameters beyond nearest neighbors in the 
dipolar spin ice model , including second neighbor ex- 
change interaction Ji- For the field- induced transition 
reported here, the role of J2 is merely to renormalize the 
value of H c slightly. 

In summary, by ac susceptibility measurements of single 
crystals of the spin ice Dy 2 Ti 2 07 under well-controlled 
magnetic fields, we finally clarified the origin of the unre- 
solved anomaly reported in the polycrystals. We clarified 
its origin as a field-induced magnetic transition governed 
by the long-range dipolar interaction beyond the nearest 
neighbors, as predicted by the recent theory by Ruff et 
al. 
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